program test
    use common_functions
    implicit none
    integer, parameter :: p=3
    double precision :: xn(p+1),w(p+1)
    call cf_lglnodes(p,xn,w)
    write(*,*) xn
    print*, w
    print*, 'xxx'


end program test